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Abstract 

The DGLAP, BFKL, modified DGLAP and modified BFKL equations are con- 
structed in a unified partonic framework. The antishadowing effect in the recom- 
bination process is emphasized, which leads to two different small x behaviors of 
gluon distribution. In the meantime, the BFKL equation and modified DGLAP 
equation are viewed as the corrections of the initial gluon correlations to the evolu- 
tion dynamics at twist-2 and twist-4, respectively. A partonic explanation of Regge 
theory for the BFKL dynamics and the relation of the BFKL dynamics with the 
helicity configuration of the initial gluons are presented. 
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1 Introduction 



The evolution equations for parton distributions play an important role in the QCD 
studies of high-energy processes. A standard evolution equation is the DGLAP equation 
(by Dokshitzer, Gribov, Lipatov, Altarelli and Parisi in [1]), which sums large logarithmic 
Q'^ corrections {LLi^Q"^)) to the integrated parton distributions. The DGLAP equation 
can be obtained by using the renormalization group or graphic techniques. The DGLAP 
equation successfully predicts the Q^-dependence of parton distributions in a broad kine- 
matic range. However, the DGLAP equation is unsatisfactory at small x, where the 
resummation of the leading logarithm of x {LL{l/x)) becomes important and it leads to 
the BFKL equation (by Bahtsky, Fadin, Kuraev and Lipatov in [2]). The BFKL equation 
for the unintegrated gluon distribution is originally derived by using the Reggeization of 
gluons. Both the BFKL and DGLAP equations predict a rapid increase of the parton 
densities at small x due to parton splitting, and the unitarity limit is violated. Therefore, 
the corrections of the higher order QCD, which suppress or shadow the growth of parton 
densities, become a focus of intensive study in recent years. Several nonlinear evolution 
equations are proposed at small x. Some of them are the GLR-MQ equation (by Gribov, 
Levin and Ryskin in [3] and by Mueller and Qiu in [4]), modified DGLAP equation (by 
Zhu and Ruan in [5]), JIMWLK equation (by Jahhan-Marian, lancu, McLerran, Weigert, 
Leonidov and Kovner in [6]), Balitsky-Kovchegov equation [7], and various versions of 
the evolution equations based on the color dipole picture [8]. Unfortunately, we still 
lack a unified picture to clearly understand the relation between the DGLAP and BFKL 
equations and the unitarity corrections on them. 

In this paper we try to present a unified partonic framework for understanding some 
of QCD evolution equations at small x. Our idea is straightforward. We begin with an 
elementary amphtude Fig. la of deep inelastic scattering (DIS) in the parton scattering 



picture, where the dashed hne imphes a virtual current probing gluon. This amphtude, 
together with its conjugate amphtude, constructs the DGLAP equation for gluons. The 
correlations among initial partons are neglected in the derivation of the DGLAP equation. 
Obviously, this assumption is invalid in the higher density region of partons, where the 
parton wave functions begin to spatially overlap. Therefore, the leading corrections of the 
correlating initial gluons to the elementary amplitude at small x should be considered. 
To this end, we add the possible initial gluons to Fig. la step by step. The resulting 
three sets of amplitudes are listed in Fig. Ib-ld. It is interesting that these amplitudes 
produce the BFKL equation, modified DGLAP equation and modified BFKL equation, 
respectively. 

A key problem in the restoring unitarity is the origin of the negative corrections. In the 
viewpoint of elementary QCD interaction, the suppression to the gluon splitting is naively 
from its inverse process-the gluon fusion (or the recombination process, i.e., a gluon split- 
ting combines with the gluon fusion) . DIS structure functions are the imaginary parts of 
the amplitudes for the forward Compton scattering of the target with a probe. According 
to the time-ordered perturbative expansion of the statement of unitarity of the S-matrix, 
the structure functions are associated with the sum of cut diagrams. These different 
cut graphs represent various possible sub-partonic processes due to the unitarity of the 
perturbative S-matrix. Therefore, the sum of cut graphs in the recombination process is 
necessary not only for infrared safety, but also for restoring unitarity. As Gribov, Levin 
and Ryskin pointed out that the 2 — > 2-2 — > 2 cut-diagram of the recombination process 
always contributes the positive correction and it can not lead to the negative effect in 
parton evolution equation [3] . The negative screening effect in the recombination process 
can occur in the interferant or virtual cut-diagrams of the recombination amplitudes [5] . 
For computing the contributions from the interference processes, work [3] used the AGK 



cutting rule [9] . The result is the GLR-MQ equation. The GLR-MQ equation was gener- 
alized to include the multi-gluon recombination process in the Glauber-Mueller formula 
[8]. 

In the next advance, one of us (zhu) disputed in [5] the above mentioned applications 
of the AGK cutting rule in the GLR-MQ equation. The TOPT-cutting rule based on 
the time ordered perturbation theory (TOPT) replaces the AGK cutting rule to expose 
the relations among various cut diagrams in a recombination process [5]. Thus, we can 
completely compute the contributions of the gluon recombination process. Following it, 
the modified DGLAP equation was proposed in [5] . A remarkable property of this equa- 
tion is that the positive antishadowing and negative shadowing in the nonlinear evolution 
equation are naturally separated. The recombination among the gluon at Bjorken vari- 
able y > X will reduce the density of the gluon at x. In the meantime, a gluon in the 
region [x/2,x] can generate a gluon at x if the former combines with its companion. In 
consequence, the modifications in the gluon distribution at x could come from either the 
screening (shadowing) effect of the fusion of two gluons at y > a;, or the enhancement 
(antishadowing) effect of the fusion of two gluons at y > x/2 and y < x. As a result, 
the corrections arising from gluon recombination at small x depends not only on the size 
of gluon density at this value of x, but also on the sharp of gluon density in the region 
[x/2,x]. Thus, the shadowing effect in the evolution process will be weakened by the 
antishadowing effect if the gluon distribution has a steeper form [5,10]. Although the 
gluon recombination re-distributes the parton densities in the above mentioned processes, 
it always keeps the momentum conservation. We note that the GLR-MQ equation also 
computes the positive contribution of the 2 — > 2-2 — > 2 cut diagram, the momentum con- 
servation of the recombination process is violated due to the mistakes in the application 
of the AGK cutting rule, where the contributions of all cut diagrams only differ by a 



numeric weight. 

We shall extend this work along the above mentioned researches. For this purpose, 
we give a short review about the applications of the TOPT in the factorization scheme. 
Following it the QCD evolution equations hsted in Fig. 1 are derived in a traditional 
Bjorken frame. We find that the antishadowing effect, which is a direct result of the 
momentum conservation, is non-negligible in the corrections of the gluon recombination 
to the DGLAP and BFKL equations. In particular, we predict that the gluon distribution 
asymptotically (abruptly) approaches to the saturation limit below (above) a moderate 
Q^-scale due to the antishadowing effect. 

It is interesting that, in this framework, the BFKL equation can be derived as a 
correction of the correlations of initial gluons at the twist-2, and the relation of the BFKL 
dynamics with the helicity configuration of the initial gluons becomes transparent. Then 
the partonic interpretation of Regge theory for the BFKL dynamics is presented. We also 
show that the similarity between the modified DGLAP equation and the BFKL equation 
in their ladder diagrams. We compare the modified BFKL equation with the color dipole 
model of the Balitsky-Kovchegov equation. We find that they have a similar physical 
picture but they are really different evolution dynamics. As is well known, the QCD 
evolution dynamics for parton densities is extensively applied in high energy physics. The 
discussions in this work present a new way to reveal some unknown properties, which can 
be regarded as a supplement to our knowledge about evolution equations. 

The paper is organized as follows. In Section 2 we discuss the separation of the 
evolution kernels from the TOPT diagrams using the semi-classical Weizsacker- Williams 
(VF— 1^) approximation [11]. In Sections 3 we discuss the connection between the DGLAP 
equation and BFKL equation. In Section 4 the modified DGLAP equation is generalized 
to include the BFKL equation. The discussions and concluding remarks are given in 



Section 5. 



2 Factorization 



The partonic picture of a high energy process is frame- and gauge-dependent. The 
appropriate choices of the coordinate frame and gauge are helpful for extracting the 
evolution kernels using a simple partonic picture. A partonic scattering model takes 
Bjorken frame, i.e., the momentum of the nucleon is 



while the virtual probe has almost zero-energy qo — {2Mv-\-Q'^)/AP and zero- longitudinal 
momentum = —{2Mu — Q'^)/AP as P — > oo, so that the momentum of the probe is 
mainly transverse to the nucleon direction. In the mean time, we work in the physical 
axial gauge: n ■ A — with n being the light cone vector 



The evolution kernel in a QCD evolution equation is a part of more complicated 
scattering diagram. In general, the correlations among the initial partons make these 
partons off mass-shell. One of the most important hypotheses in the derivation of the 
DGLAP equation is that all correlations among initial partons are negligible during the 
interaction. In this case, the interaction of a virtual probe with the nucleon can be 
factorized as the nonperturbative parton distribution and hard probe-parton matrix in 
the coUinear factorization scheme. 

On the other hand, at the higher density range of partons, these correlations among the 
initial partons can no longer be neglected (Fig. 2). In this case, the transverse momenta 



P'^ = (Po,Z,P3) = (P,Q,^), 



(2.1) 




(2.2) 



and 




(2.3) 



of the initial partons are non-zero and these partons are generally off mass-shell, therefore, 
the /cT-factorization scheme [12] is necessary. 

However, in this work we use the W—W approximation rather than the /cr-factorization 
scheme. The reason is that the W — W approximation allows us to extract the evolution 
kernels and to keep all initial and final partons of the evolution kernels on their mass- 
shell. For illustrating this idea, we consider a sub-process as shown in Fig. 3 , where the 
momenta of three massless particles are parameterized as 



Pi = (p,o,p), 



P2 = {X2P+- ,k,X2p), 

2X2P 

and 



P3^ipl,-k,X3p), (2.4) 

where p — > oo. The particle 3 is off its mass-shell. We decompose this propagator into a 
forward and a backward propagator according to the TOPT. The corresponding momenta 
are on mass-shell 



^2 

Pf = {xsp + -k, X3p), 

2X3P 



and 



^2 

Pb = (-X3P - TT^, -k, xsp), (2.5) 

2X3P 



respectively. Note that Eqs. (2.4) and (2.5) require 



xp:^\k\, (2.6) 

which is consistent with the conditions of high energy- and small angle-scattering in the 
W — W approximation. One can find that the forward propagator contributes the power- 
factor l/kj_, since 



while the backward propagator gives the contributions ~ l/p^- If the dominant contribu- 
tions in a process are from the terms with the power of transverse momentum, one can 
find that the contributions from the backward propagator are suppressed at p — > oo, and 
the dominant contribution is from the forward component. In this case, we can break the 
forward propagator since it propagates on mass-shell. 

The wave function of a parton in the Bjorken frame occupies a space with transverse 
size ^ 1/k and longitudinal size ~ 1/xp in the nucleon. The partons are regarded as inde- 
pendent if the partons are dilute and the parton distributions obey the DGLAP equation 
(Fig. 4a). However, this evolution equation predicts a rapid rise of gluon multiplicities 
inside the nucleon at small x. At sufficiently large gluon densities, the wave functions 
of the gluons overlap with each other. We call such a correlating gluon cluster with a 
transverse scale a;„(, as the cold spot, which phenomenologically describes the correlation 
among initial partons (the dashed circles in Fig. 4b). In this work we mainly consider 
the cold spots consisted of two gluons, which dominates the leading small x processes. 

The elementary correlation among two initial partons in a cold spot is shown in Fig. 5, 
where the dark circle implies all possible QCD-channels. Let us consider a configuration 
of two correlating partons that form a cold spot inside a nucleon. We assume that the 
longitudinal momentum of the initial parton is much larger than its transverse momentum 



in the Bjorken frame even at small x. Thus we can take the W — W approximation and 
extract the two parton correlation function from the TOPT-diagram as 



2Ef 



1 



- - Epb 



2E„ 



2E, 



Pb 
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(2.8) 



Using Fourier transform with respect to the transverse momentum transfer, these 
partons localize in a range with the size on the impact space. The unintegrated 
distribution /(x^j^; Xa, Xb) can be understood as the probability of finding two partons with 
transverse distance x^i, and the longitudinal momentum fractions Xa and Xb, respectively. 

The probabihty of forming a cold spot is proportional to the "size" of gluon. In 
particular, we have a straightforward bound condition 



fi^b = 0> ^a, Xb) = 0, 



(2.9) 



which implies that the overlap probability of two sizeZess partons in the impact space 
vanishes. 

In the impact space, the integrated gluon distribution reads 



G{xl2,x) = xg{x^2,x) 



^o2 ^ab 



/(^6,a;), 



(2.10) 



or 



2 dxg{x^2,x) 

OXq2 



—q2 —ab 



(2.11) 



while in the momentum space, they are 



G{Q',x)=xg{Q^x) 



(2.12) 



or 



fit ,x)=Q — Q2=fe2 , (2.13) 

where Rn is an effective radius of the nucleon. Note that the momentum Q of the probe 
is mainly transverse to the nucleon direction. 

One can imagine that the cold spots fill the whole transverse plane of a nucleon at 
a higher density scale Ug, which corresponds to the saturation of the gluon distributions 
(Fig. 4c). The saturation is a comphcated state, where the long- and short-distance 
correlations among partons coexist in a nucleon. In this work, we focus on the region far 
from the saturation limit. 

In the next step, we discuss the separation of the virtual probe-vertex from the DIS 
amplitude. In principle, in- and out-probe lines can be attached to the left- and right- 
hand sides of the cut line in all possible ways. However, in the Bjorken frame, one 
only need to pick up the diagrams with two photon lines following the cut line at the 
leading logarithmic (Q^) approximation {LL{Q'^)A), since these diagrams contribute the 
dominant terms ~ dh^/^ [5]. 

Under the same (i.e., LL{Q'^)A) condition, the contributions of the backward compo- 
nents in the two propagators connecting with the probe are suppressed. Thus, we can 
take the W — W approximation to factorize the virtual probe- vertex from the cut diagram 
and extract the evolution kernel in TOPT. However, the BFKL dynamics works in the 
region beyond the LL{Q^)A. Fortunately, we could show that the dominant contributions 
still come from the terms with the power of the transverse momentum, i.e., we pick up 
the dominant terms ~ d^/k"^. Therefore, the W — W approximation is still feasible in 
the separation of the probe-vertex for the BFKL dynamics. 



3 Evolution equations at twist-2 

Using the W — W approximation mentioned in Sec. 2, the DIS cross section at leading 
twist is generally factorized into the following convolution form 



da{probe*p — >• k'X) 



fikl,xi)0}C 




(3.1) 



with the perturbative evolution kernel /C, the pro6e*-parton cross section da{probe*k2 — > 
k') and the nonperturbative unintegrated parton distribution function (s) /. For simplicity, 
we take the fixed ag in this work. Wc emphasize that the factorization between IC and 
da{probe*k2 — > k') is valid because the propagators linking these two parts are on mass- 
shall in the W — W approximation, and hence we can define Af{h^;x2)- According to 
the scale-invariant parton picture of the renormalization group [13], we regard Af{^, X2) 
as the increment of the distribution f{^,Xi) when it evolves from {^,Xi) to {^,X2)- 
Thus, the connection between the unintegrated gluon distribution functions f{k^, Xi) and 
f{kl,X2) due to Eq. (3.1) is 



where the factorization scales are integrated, thus the observed cross section da{probe*p — > 
k'X) is independent of the factorization scheme. 

In the case of the evolution along the transverse momentum, we differentiate Eq. (3.2) 



f{klx2)^f{klxi)+Af{klx2) 



f{klxi) + 




(3.2) 



with respect to and get 



df{kl,X2) 



/•I dXi 1 (kl X2 _\ , 2 ^ 



^1 



+ L (3.3) 

Unfortunately, this equation is unsolvable in a resummation form. However, at the 
LL(^)A, the evolution kernel JC in Eq. (3.3) only is the function of the longitudinal 
variables and the second term in the right-hand side of Eq. (3.3) vanishes. In this case, 
using Eq. (2.12) we have 



AG(g^x,)^ /f ^Af{klx2) 



dB^ rkl dkl d^i ^ _ 



u2 t,2 / 2 J,2 

!k2min J±.2 -Imin 1^ 

^ dk^ dxi X2 



h r dXi (X2 _ \ , 2 ^ 
I Jx2 Xi \Xi J 



dk^ f'dXiX2 X2 _\ 2 ^ (o .^ 

L T2- / ^DGLAP — , as G[k2, Xi , (3.4) 

4Lin &2 -^^2 Xi Xi \Xi J 



and 



G{Q\ X2) = G'(^^, xi) + AG'(g^ X2), (3.5) 
which lead to the DGLAP equation for the gluon distribution 

^dg{Q^,X2) dxi (^^ -\ ^^2 x /q 

^ = / f^DGLAP [—,as] g{Q , Xl). (3.6) 

OL^ Jx2 Xl \Xi / 

We could extract the evolution kernel JCdglap in the coUinear factorization scheme [14] 
if the transverse momenta of the initial partons are neglected. 

It is interesting that the DGLAP equation Eq. (3.6) also can be derived from the 
renormalization group, where the gluonic structure function is factorized as 



F(g^ X,) = /' ^, as{A g(^', x,). (3.7) 
Jx2 Xi Xi J 

In this equation /i is the unintegrated renormahzation scale as well as the factorization 
scale, C is the perturbative coefficient function. Using 

dF(Q'^,X2) , , 

' = 0, (3.8) 

one can get [15] 



/■I dxi 

'x2 Xl 



din fjL^ 



g[li , + C -, — , a, /X — — 

\li Xl J din/i^ 



0, 



(3.9) 



at the lowest order of ^^(/x^), it becomes 



Jx2 Xl yn Xl J dm/i^ 



-L 



^dxidC^'^ 



X2 Xl 



din. 



g{fJ'^,Xi). 



(3.10) 



Using 



C(o) = 5 (i - ^) , (3.11) 

and 

C^'^^ICGDLAp(-)ln%, (3.12) 

one can get the same results as Eq. (3.6). 

One difference between the coefficient function C^^-* in Eq. (3.10) and the evolution 
kernel JCdglap in Eq. (3.6) is that two propagators connected with the probe in the co- 
efficient function C^^^ are off mass-shell, but not on mass-shell hke in Eq. (3.1). However, 



as we have proved that the evolution kernels from two methods are really equivalent in the 
infrared safe theory even at twist-4. The reason is that the contributions of the backward 
components in these two propagators always are suppressed [5]. 

Considering that the scale- invariant parton picture [13] is a phenomenological inter- 
pretation of the renormalization group in the derivation of the QCD evolution equation, 
the above mentioned two methods are really equivalent. 

We now construct the evolution equations along the longitudinal variable. As we shall 
show that the transverse position of the targeted gluon is uncertain in this case, we need 
to change the integrated phase space in Eq. (3.2), i.e., 

r dkl r dxi r dkl r dx2 kl X2 , . 

J kl J xi J kl J X2 tlxi 

Thus, we have 



f{klx,) = f{klx2)-Af{klx,) 




(3.14) 



where we omit the integral limit for the transverse momentum since they are non-ordering. 
At the LL{l/x)A, K,' is only the function of the transverse momenta. We differentiate 
Eq. (3.14) with respect to xi and obtain 




This is the real part of the BFKL equation. Where '— ' imphes that the evolution along 
the direction of decreasing — dx. 

Now let us calculate the evolution kernel. We take the square of the total amplitude 
in Fig. lb and separate out the probe vertex using the W — W approximation [11], one 
can get the four TOPT-diagrams Fig. 6, where the evolution kernels are constructed by 
following amplitudes: 



^BFKL — ^BFKLl + ^BFKL2i 




and 




(3.16) 



The momenta of the partons are parameterized as 



Pa = {xiP + 



,k + l„xiP), 



(3.17) 



2xiP 



k = {X2P + 



2X2P 



,k, X2P), 



(3.18) 



la ^ {{Xl - X2)P + 



2{xi-X2)P 




l_a,{xi- X2)P), 



(3.19) 



Pb = {xiP + 



,k + lj„xiP), 



(3.20) 



2xiP 



and 



;2 

h = {{xi - X2)P + ^, .^ ,h, (xi - X2)P). (3.21) 



One of the matrix is 



MbFKLI = i9r''''[gafi{Pa + k)^ + gp^{-k + la)a + g^a{-la-Pa)l3]^a{Pa)^l3{k)e^{la), (3.22) 

where the polarization vectors are 



6(p„)^(0,6,- -"^- + -"\ (3.23) 

and 



eW = (0,e,-|-|), (3.24) 



Taking the leading logarithmic (1/x) approximation, i.e., assuming that X2 -C Xi, one 
can get the total amplitude 

ABFKL{k^,,k^^n,xuX2) = tgr'^2M[^ + ^-^^ (3.26) 

V ^2 Kup ^down 

where we use different indices up and down to distinguish between the different paths of 
the probing position along the time order on the impact space. 

The impact parameter in the evolution equations is first introduced by A. Mueller in 
his color dipole approaches [8]. We imitate his method but work in the scattering model. 
After Fourier transforming with respect to the transverse momentum transfer, we obtain 
the representation of the amplitude in the impact parameter space as follows 



ABFKL\XaQ,Xs;}h,Xl,^2) — J (27rj^ ^BFKL\^upi ISulowm ^2)^ 



= igr'^2J :ii - =1^] . e. (3.27) 

V ^2 XaO ^bO 

Using this amplitude we construct the evolution kernel in the BFKL equation 



:< 3 '>colar 2~2 2"*^ — ' (3.28) 
TI" ^0^0 ^2 



where the position Xq in Eq. (3.28) is running along the transverse coordinator due to 
uncertainty relation although k^^ identifies with ^own forward deep inelastic scat- 

tering. Therefore, the measurement of the unintegrated parton distributions in the BFKL 
dynamics is non-local and we use the dark box to indicate the above mentioned measure- 
ment. On the other hand, we can not regard the fraction x of the longitudinal momentum 
of the parton in the BFKL equation as the Bjorken variable since the transverse scale ]^ 
of this parton is irrelevant to the probe scale Q^. 
Using Eqs. (3.15) and (3.28), we get 



df{xl^,x) 



dx 

^OLg f j2 cib J? / 2 \ 

, 2 / d'x,^[f{xl,,x) + f{xl,x)]. (3.29) 



3a 



Note that the distribution /(xq, x) of a gluon at inside a cold spot, which is scaled by 
a fixed Xai, can be represented by either f{x^Q, x) (if taking the coordinate relative to Xg), 
or f{^j„x) (if taking the coordinate relative to Xj,). The right-hand side of Eq. (3.29) 



is illustrated in Fig. 7b, where the un-connecting gluons in Fig. 7 means all possible 
connections as shown in Fig. 6 and in the last diagram the probe is included. 

The equation (3.29) is the real part of the BFKL equation in the impact space. A 
complete BFKL equation should include the contributions from the virtual processes, 
which are necessary for the infrared safety of the evolution dynamics. Using the TOPT- 
cutting rule (see Appendix), one can prove that the virtual diagrams in Fig 8 (where we 
omit the conjugation diagrams) contribute the same evolution kernel as the real kernel but 
differ by a factor —1/2. In the mean time, the probed distribution f{x^Q, x) is replaced by 
/(x^j, x) because of the change of the cut hne. Thus, we can directly write the complete 
equation as 



—X- 



df {xlf„x) 
dx 



27r2 



f (^aO ) ^) + / (^60' x) f (^6' x) 



(3.30) 



This form of the BFKL equation was first derived by Mueller in the color dipole model 
[8]. 

To make regularity of the evolution kernel explicit, we use (for simplicity we replace 

^ao ^ Pi and x^Q ps) 



The evolution equation becomes 



Pi P2 P1P2 



(3.31) 



X 



dfjp, x) 
dx 



^ 3a, r 

TT Jo 



dpi 



\Pi - P\ 



—f{Pi,x) - f{p,x) 
Pi 



3as dpi 

TT Jp pi 



(3.32) 



The first integral of Eq. (3.32) is regular both for pi — > and pi —>■ p. Defining pi = up 
when pi < p and pi — p/u when pi > p, one can write Eq. (3.32) as 



dx 

= — /' T^ifiup, x)/u + f(p/u, x) - 2/(p, x)]. (3.33) 

TT Jo 1 — M 

The solution at small x is power behaved 

f{di.,x)r^^/^x-\ (3.34) 

where 

A = 4^5 In 2. (3.35) 

Now let us discuss the relation of the BFKL dynamics with the DGLAP evolution 
equation. The correlations among the initial gluons disappear in the dilute parton system. 
In this case the contributions of the interference diagrams Figs. 6c and 6d are negligible. 
Thus, Eq. (3.28) reduces to sum two independent sphtting functions, each of which is the 
splitting function in the DGLAP equation in the impact space at the small x limit, 

dx^dxi sr^ I i I \\i d^la 



^0 ^1 pol )^^la 



dxi dx^ 
27r X2 xlo 



■aO 



(3.36) 



Because the interference amplitudes disappear in the DGLAP dynamics, the transverse 
position of the probing parton is irrelevant to the measurement, a simple probe (dashed 
line in Fig. 7) can be used. On the other hand, a single transverse scale x^q in the DGLAP 
dynamics allows us to take = Thus we can define the Bjorken variable xb and set 
X2 — Xb in Eq. (3.36). Returning to the momentum space and using Eq. (3.6), one can 



get the real part of the DGLAP equation at the double leading logarithmic (i.e., InQ^ 
and ln(l/x)) approximation {DLL A) 



dg{Q^,XB) 
91ng2 



/ —}CDGLAp{z,as)g{Q'^,xi)--g{Q'^,XB)[ dzK,DGLAp{z,as) 





(3.37) 



where the contributions of the virtual diagrams disappear at small x since they company 
with the singularities at a; — > 1 (see appendix A). 

The above derivations also reveal a relation between two evolution equations from a 
new view point: the evolution kernel of the DGLAP equation is the un-interferant part 
of the BFKL kernel. Or reversely, the BFKL equation can be regarded as the corrections 
from the initial parton correlations to the DGLAP equation on the twist-2 level. 



4 Evolution equations at twist-4 

We consider the evolution kernel based on Fig. Id and call it as the modified BFKL 
equation. Note that two pairs of initial gluons, which are hidden in the correlation func- 
tion, for example in Fig. 9a, should be indicated as Fig. 9b. A set of cut diagrams based 
on the Fig. Id are listed in Fig. 10, where the probe vertices have been separated out 
using the W — W approximation. Similar to the derivation of Eq. (3.28), we write the 
evolution kernel in this equation as 

ICmd-bfkl = X] ^md-bfklA\jd-bfkl ^2'jr^)2Ei ' ^^''^^ 

The amplitudes 



■^MD-BFKL — ^MD-BFKLl + ■^MD-BFKL2, (4-2) 

where 



and 



^MD-BFKLl 



I Ek 1 1 

Epa + Ep^ 2Ek Ek + Ei^ - Ep^ - Ep^ 



MD-BFKLl, 



(4.3) 



^MD-BFKL2 



Eu 



-M. 



MD-BFKL2- 



i Ep^ + Ep^ 2Ek Ek + Ei^ - Ep^ - Ep^ 
For simplicity, in Fig. 10 we assume that all cold spots have a same size 



(4.4) 



^6 



X 



X 



c'd' 



and 



— Sj)'c'- 



(4.5) 



The momenta of the partons, for example, are parameterized as 



Pa = {XiP + :L-m: X^P), 



( _|_ fYi\^ 

Pb = {xiP+ ~ ~ ,k + rn,xiP), 

ZXiP 



kb = {X2P + ^^,k,X2P), 



i. = ((2.,-.,)P+^p-4-^.l,(2.,-.,)P), 



and in the t-channel 



(k + mf k^ . 
m^Pb-kb^ {{xi - X2)P + ^^^p ^^^,m-{xi -X2)P). 

The matrix in Eq. (4.3) is 



Mmd-bfkli = igr''C^''''^^tgf'^Cf''''^eMepiPbyp{la)e:{k), 
where C°'(^'^Cf"^'^ are the triple gluon vertices and polarization vectors are 



e{Pa) = {Q,e,- - 7p-0 , 



/X / e ■ ik + m) , 

6(p,)^(0,6,- - ^-p-^ ), 



e ■ k 



and 



Thus, we have 



(4.15) 



AMD-BFKL{k : kdown^ 221, 22l', Xi, X2) 



= /r-r-./l^ie -'f;"^-;- + 6 ^-ff-^;fV (4.16) 

Taking the Fourier transformation 

AMD-BFKL{s.m^ ^c, ^ab) ^cd, ^1, ^2) 

= / ^'^»^^'^--f ^^'^^ ^^^_^^^,(fe^^, ^, rd. x„ x,)e-^^.^^°+^^c,..„.oc+^ir^..+^M'.cM 

(4.17) 

we get 



/I / \ /-■ 2 £abc £dce / "^l r— WSIafe 2Ic02iafei //i io\ 

>lMD-BFi^L(^60,^cO,^a&,a;i,X2) = /""V Wt^I^^ • 6. (4.18) 

Summing all channels, we derived the evolution kernel corresponding to Fig. 10 as 



jqMD-BFKL ^2^6 ^^^2 



2 2 
^0 ^6 ^2 



1 c?a;2 2.. .,2 

.^MD-BFKL^MD-BFKL 
poi L-LOTT X2 



^MD-BFKL^*K 



2(-^) 36q;^ ^^^^6c ^^^b ^^^^^ 

32 coior TT ^^0— cO ^ab "^2 



In the case of \xjj^\ ^ can neglect the correlations between two cold spots 

with the sizes x^j, and x^,^. In consequence, the contributions of the interferant terms 
(Figs. 10c and lOd) disappear and Fig. Id backs to Fig. Ic. Thus, Eq. (4.19) reduces to 
the real part of the modified DGLAP-kernel 



j^MD-DGLAP ^^O '^^1 _ 2^ ^ ^ SGctg (i^Xfep dXj (4 20) 

—60 "^1 color TT 2^2 

Similar to Eq. (3.5) we derive the modified DGLAP equation using 



G{Ql X2) = G{Ql xi) + AG{Ql X2) 

= GiQlx,) + ^ /"^' ^^X:md-z.glap (^,as) G^'\qIx,), (4.21) 

where a power suppressed factor 1/Ql has been extracted from the evolution kernel. The 
4-gluon correlation function G^"^^ is a generalization of the gluon distribution beyond the 
leading twist. It is usually modelled as the square of the gluon distribution [3,5]. For 
example, G^^^ = CiG^ and Ci = l/{7rRjf). The evolution constructed by ICmd-dglap is 
illustrated by Fig. 7c. 

The complete modified DGLAP equation combining DGLAP dynamics at small x and 
in the momentum space was written as [5] 

dGjQ^XB) 
91ng2 

where the right-hand second term of (4.22) is positive antishadowing and comes from the 
contributions from Fig. 11a, while the third term is negative shadowing and arises from 
the contributions of the 1 — > 2-3 — > 2 cut diagrams (see Fig. lib, which is the same 
order as Fig. 11a but we haven't indicated it in Fig. Ic). Note that the shadowing and 
antishadowing terms are defined on different kinematics domains [x, 1/2] and [x/2, 1/2], 
respectively. Comparing with the GLR-MQ equation, there are several features in the 
modified DGLAP equation: (i) The momentum conservation of partons is restored in a 



complete modified DGLAP equation; (ii) Because of tlie sliadowing and antisliadowing 
effects in tfie modified DGLAP equation fiave different kinematic regions, tfie net effect 
depends not only on the local value of the gluon distribution at the observed point, but 
also on the shape of the gluon distribution when the Bjorken variable goes from x to 
x/2. In consequence, the shadowing effect in the evolution process will be obviously 
weakened by the antishadowing effect if the distribution is steeper [5]; (iii) Both the 
GLR-MQ and modified DGLAP equations assume that the fusion only occurs among 
gluons with the same value of x. However, there is a priori no reason to forbid the 
recombination of two gluons with different values x. If the recombination of gluons with 
different x is included, we find that this modification unreasonably enhances the shadowing 
effect in the GLR-MQ equation, while it does not change the predictions of the modified 
DGLAP equation. Therefore, a correct antishadowing correction is important in gluon 
recombination [10]; (iv) Finally, we emphasize that, as we have shown in [10], the existence 
in the recombination process is a general result of the momentum conservation and it is 
irrelevant to the concrete form of the evolution dynamics. 

Now let us discuss the modified BFKL equation using Eq. (4.19), which is the cor- 
rections of gluon recombination to the BFKL equation and it is illustrated in Fig. 7d. 
Unfortunately, it is difficult to write a simple form of the modified BFKL equation due to 
the complicated evolution kernel (4.19). For the sake of simplification, we take the follow- 
ing approximations (see Fig. 12): At the first step, we take the W — W approximation to 
factorize Fig. 12a to a four gluon correlation function (Fig. 12b) and the general recom- 
bination function (Fig. 12c). Comparing Eq. (3.36) with Eq. (3.28), one can find that 
the kernel JCdglap is a part of the kernel JCbfkl and they differ only by a phase space. 
Similarly, the modified DGLAP kernel is evolving along the transverse momentum in the 
modified DGLAP equation, but, as a part of the modified BFKL-kernel (Fig. 10a and 



10b), they also participate the evolution along the longitudinal momentum. Therefore, 
at the next step, we only keep the contributions of Fig. 10a in Fig. 12c, i.e., using the 
modified DGLAP kernel to replace the kernel (4.19). According to Fig. 12, we rewrite 
l^G{Q2,X2) in Eq. (4.21) in the impact space as 



AG(xL,x.) =j/-^j ^^ICmd-dglap (^,as) f^'\^,x,). (4.23) 
Thus, we have 



dA G{xlf;„X2] 
dxi 



iab -60 



= \ r —-^MD-DGLAP (^,«.) f'\x^,.X,), (4.24) 
X_gjy Jxiarxill X\ X\ \X\ J 

which gives a correction to the evolution of the unintegrated distribution along small x 
direction 



dfjx^h^X) _ 81 r(2)( 2 ^^ ^^'^sr(2)( 2 x ( a 



Combining with the BFKL equation, we obtain an approximate form of the modified 
BFKL equation 



^ df{xlj,,x) 
dx 

Oia f ,0 „ x^h o N 1 Q; 



/ d?x^Q^f{xi,,x) - ^^/(x^,,x) / £x,6- 

J X„nXi.n Z ZTT J 



(4.26) 

where the nonlinear evolution kernels are the derivative of the modified DGLAP kernels 
with respect to x; in the meantime, we assume that /^^^ = (^2/^ and C2 is a unknown 
correlation coefficient. 



Comparing with the modified DGLAP equation (4.22), the antishadowing effect in the 
modified BFKL equation (4.26) is more sensitive to the shape of the gluon distribution 
at small x. For example, in the case of a typical BFKL solution Eq. (3.34), the anti- 
shadowing effect may cancel, or even outwegh than the shadowing effect in the modified 
BFKL equation (4.26). Such faster increasing gluon distribution impels the multi-gluon 
recombination (or correlation of multi-cold spots) to participate the evolution and reach 
the saturation limit earlier. 

There is a difficulty in the numerical calculation of Eq. (4.26): the value of f{x^i„ x/2) 
is unknown when the equation evolves to small x. To avoid this problem, we take following 
approximation at the beginning modified BFKL-renge, for example, in the evolution from 
xi to X2 — xi — Ax, the increment of the unintegrated gluon distribution due to Eq. 
(4.26) is taken as 



where BFKL, AS and S indicate the contributions from the BFKL kernel, antishadwoing 
term and shadowing term in Eq. (4.26) and /efkl imphes the solution of the linear BFKL 
equation. A schematic solution of Eq. (4.26) at this approximation is illustrated in Fig. 

13. Where we take a typical BFKL input distribution as 



and = IGeV^, 

3^b ~ ^GeV ^. Comparing the coefficients of the nonlinear terms of 
Eq. (4.26) with that of Eq. (4.22), we take C2 = l/{25n)GeV~^. 

The dashed-point curve in Fig. 13 is the imaginal saturation limit. One can find that 
the antishadowing effect is un-negligible in such abrupt distribution. For comparison, we 



Af{xl„ X,) = BFKL[f{^„ xo)] + AS[f{^„ -j)] - S[f{c^„ xo)] 




(4.27) 




(4.28) 



draw a corresponding BFKL-solution (dashed curve) and a solution of Eq. (4.26) without 
the antishadowing corrections (point curve) in Fig. 13. 

The schematic kinematic ranges of four relating evolution equations are shown in Fig. 
14. Although both the BFKL equation and modified DGLAP equation are regarded as the 
corrections of the correlations among initial partons to the DGLAP dynamics, the BFKL 
equation dominates the higher Q^-range since the modified DGALP equation is order of a"^ 
and has a power suppression factor. However, the modified DGLAP equation will replace 
the BFKL equation below a moderate Q^-scale because the 1/x factor in KLmd-dglap 
(see Eq. (4.20)). It is interesting that because the gluon distribution becomes fiatter in 
the modified DGLAP-region in Fig. 14, the antishadowing contributions in Eq. (4.26) can 
be neglected when the evolution enter into the modified BFKL-region. Thus, we predict 
that the gluon distribution asymptotically or abruptly approach to the saturation limit 
below or above a moderate Q^-scale. 



5 Discussions 

At first, we point out that our nonlinear equation (4.26) has a similar picture of the 
Bahtsky-Kovchegov equation [7], but they have difTerent interpretations. The Balitsky- 
Kovchegov equation is one of nonlinear evolution equation incorporating with the screen- 
ing corrections to the BFKL equation. Besides, this equation is regarded as an approxima- 
tion form of the JIMWLK equation at the lowest order. According to Ref . [8] , assuming 
that the scattering matrix 5'(x„^,x), which is the part of dashed box in Fig. 7b but use 
quark- antiquark dipole to replace the gluonic dipole, obeys the BFKL equation 



dS{xlt„x) 



dx 

- / (fxQ— ^[S{a^Q,x) + S{a^Q,x) — S{a^jj,x)]. (5-1; 



27r2 

Then assuming that 

(1) The gluon decay inside the dipole leads to the dipole splitting at large Nc, i.e., 

S{x^o, x) + S{2^, x) — > S{x^Q, x)S{xIq, x), (5-2) 

in Eq. (5.1); 

(2) Defining the scattering amplitude 

S =1-T, (5.3) 
one can simply obtain the Balitsky-Kovchegov equation 



/T{xl„x) 
dx 

/ d\,^^[T{^^, x) + T{^, x) - T{^„ x) - T{xl„ x)T{xlo, x)], (5.4) 



where T is relative to the gluon distribution [7]. Figure 15 is the one-step evolution 
containing the 4 — > 2 recombination amplitude. 

Comparing Fig. 15 with Fig. 7d, one can find that the Balitsky-Kovchegov equation 
and the modified BFKL equation have a similar physical picture. However, they are 
really different nonlinear evolution dynamics. In fact, the replacements Eq. (5.2) and 
(5.3) result the negative term TT in Eq. (5.4). This term is graphically regarded as 
the simultaneous scattering of two dipole. Such consideration of the shadowing effect is 
relevant neither to the AGK- nor to TOPT-cutting rules, but it is similar to a shadowing 
mechanism in the multi-scattering process (for example, the Glauber-type scattering ). 

Different from the recombination process, the antishadowing effect is absent in the 
Glauber-type scattering since the shadowed (or absorbed) momenta by target are un- 
measured [17]. Besides, the linear parts and nonlinear part in the Balitsky-Kovchegov 
equation share the same BFKL-kernel due to Eqs. (5.2) and (5.3). On the other hand, 
the nonlinear evolution kernel in the modified BFKL equation (4.26) is different from 
the linear BFKL-kernel since the former including the recombination process. Thus, we 
present two different shadowing mechanisms: the QCD gluon recombination process and 
the absorptive effect in the multi-scattering process. 

We have shown a new approach to derive a set of consistent evolution equations in 
a partonic framework. We present the connections among the DGLAP equation (3.37), 
BFKL equation (3.30), modified DGLAP equation (4.22) and modified BFKL equation 
(4.26) in a clear physical picture (Fig. 1). We also show that these evolution equations 
are uniformly derived from the factorizability of DIS cross sections. The four evolution 
equations mentioned above have similar structures: The positive part in an evolution 
equation is completely separated from the negative part. The positive part is the contri- 
butions from the real diagrams, while the negative part comes from the virtual diagrams 



in Eqs. (3.30) and (A-9), or from the interference diagrams in Eqs. (4.22) and (4.26), 
respectively. 

The BFKL equation is traditionally understood as an improvement of the DGLAP 
equation at small x, where the summation of the leading logarithm of x {LL{l/x)) terms 
is important. The derivation of the BFKL equation in this work shows a new relation 
with the modified DGLAP equation. We re-draw the cut diagrams Figs. 16a- 16d in our 
derivation of the BFKL equation as the ladder diagrams in Figs. 16e-16h. On the other 
hand, a typical diagram in Regge approach of the BFKL equation is shown in Fig. 17, 
where the vertical (dashed) lines are reggeized gluons, which can be explained as the gluon 
ladders in perturbative QCD. Thus, Figs. 16e-16h give a possible partonic explanation of 
the Regge theory for the BFKL equation. 

As is well known, the BFKL equation could also be derived in the color dipole model 
[8]. Now let us compare the evolution equations in our partonic model with that in 
the color dipole model. We point out that these two approaches have different partonic 
configurations. The dipole is constructed by two heavy quarks, which represent a whole 
nucleon in the color dipole model, while the cold spot in our approach is generally the 
multi-parton cluster (as one of the constituents of the nucleon). In the concrete, we 
compare Fig. 18a in our partonic scattering model with the color dipole picture Figs. 18b 
or 18c, where the transverse coordinators of the dipole are assumed to be freezed during 
interaction. Note that the dashed lines in Fig. 18a are integrated out as the final state in 
the inclusive DIS processes, therefore, some of them are omitted in Fig. 1. On the other 
hand, the splitting kernel in the color dipole model is constructed by the gluon-quark 
vertex, while it is the three gluon vertex in our partonic scattering approach. The latter 
dominates the small x behavior since it (~ 1/2:2). 

It is interesting that a cut diagram Fig. 19a of the modified DGLAP equation is 



equivalent to the ladder diagram Fig. 19b. Comparing Fig. 19b with Figs. 17e-17h, we 
find that the BFKL and modified DGLAP equations have the same (two ladders) structure 
in the initial gluon configuration, but the probe partially couples with two ladders in the 
BFKL equation at twist-2, while it fully couples with two ladders though the gluon fusion 
mechanics in the modified DGLAP equation at twist-4. 

We now discuss the relation of the BFKL dynamics with the hehcity configuration of 
the initial gluons. As we have pointed out in [5], the polarized form of two correlating 
initial partons is relevant to the structure of the spin-independent evolution kernel. For 
ease of representation, we draw all possible helicity amphtudes in Fig. 20. Obviously, the 
processes illustrated Figs. 20(al)-20(a4) should be inhibited due to all initial and final 
states in Fig. 20 having fixed helicities in the W — W approximation. Thus, the BFKL 
kernel (3.22) is a result, where the contributions from Figs. 20(bl)-20(b4) being excluded, 
i.e., we should assume that two gluons in a cold spot possess the same hehcity. On the 
other hand, the BFKL kernel reduces to the DGLAP kernel if these two correlating gluons 
have the opposite helicities, or the BFKL kernel becomes 



dx^ dx2 

BFKL— 2 

^0 ^2 



as 

"■^ >J color 9 



+ 



2 dX2 

d Xn , (5.5) 

X2 



if two gluons in a cold spot can take any helicities. 

In summary, we discussed the corrections of the initial gluon correlations to the evolu- 
tion dynamics in a unified framework. We presented the following results: (1) a possible 
connection among a set of evolution equations at small x; (2) the antishadowing effect, 
which is a direct result of the momentum conservation in elementary QCD process, is un- 
negligible in the corrections of the gluon recombination to either the DGLAP or BFKL 
equation; (3) the gluon distribution asymptotically (abruptly) approaches to the satura- 



tion limit in lower (higher) regions is predicted; (4) the BFKL equation and modified 
DGLAP equation are viewed as the corrections of the initial gluon correlations to the 
evolution dynamics at twist-2 and twist-4, respectively; (5) a partonic interpretation of 
Regge theory for the BFKL dynamics and the relation of the BFKL dynamics with the 
helicity configuration of the initial gluons are presented. 
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Appendix: 

For illustrating the applications of the TOPT-cutting rule, which are proposed in [5] , 
we write a complete derivation of the DGLAP equation for gluon. The contributions of 
the real diagrams (Fig. 21a) to the hard part of Eq. (3.1) are 
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1 1 1 1 1 

Mi^ki' 

2Ei 2Ek El — Ek — Eii 2Ek Ei — E^ — Ei 



-MUv 



dH' 



{2ttY2Ev 



1 d^k' 

X— \Mprobe*k^k'\'^{'2'^)'^S{Pprobe* + k- k')- 



4E, 



probe 



{2t:Y2Eu' 



dp 

= ^DGLAp{xi,X2, xs, as)dx:i—^5{x2 - X B)dx2C{prohe*k k'), 
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where the factor l/{2Ei) is extracted from the definition ofthe gluon distribution g{Q'^,XB), 
and 
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Mprobe*k^k'\'^{^^yS{Pprobe* + k - k') 



d^k' 
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is the contributions from da{probe*k — > k'). Thus, 
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where we inset 5{xi — X2 — X3)dxi. We obtain 
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On the other hand, the contributions from one of the virtual diagrams, for say, Fi 
(21b) are 
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is the contributions from da{probe*l k'). Therefore, we have 
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^ -^g{Q^,XB) J dz)CDGLAp{z,as) . {A 
The complete evolution equation for the gluons is 
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Note that the factor 1/2 in Eq. (A-8) is the considerations of the symmetry under ex- 
change of two internal gluons in the virtual diagrams (see Ref. [5]). The expressions 
of the evolution equation in the TOPT form, without the calculations of the matrixes, 
show that the real and virtual diagrams contribute the same evolution kernel but with 
the different factors. This simply form of the equation consists with a more complicate 
derivation of the DGLAP equation using the covariant perturbation theory in [19]. In 
fact, from Eqs. (A-9) and (A-10) we have 
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Note that dzz/ {1 ~ z) = Jq dz{l — z)/z, we obtain 
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This equation is equivalent to the following traditional for of the DGLAP equation 
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Equation (A-14) becomes 
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Both Eqs. (A-12) and (A-16) at small x predict Eq. (3.37). One can find that the 
infrared divergence of the real diagram at the end point {xb 1) is exactly cancelled by 
the the contributions of the virtual diagram. Therefore, the leading terms of the virtual 
diagram in Eq. (3.37) disappear at the small x region. It is different from Eq. (3.37), the 
contributions of the virtual diagram exist in Eq. (3.30), which regularize the singularities 
of the real diagram on the transverse space. The above results also can be apphed in the 
BFKL dynamics. A difference in this case is that the absence of the conditions 6{xi — xb) 
and 5{x2 — Xb) in Eqs. (A-1) and (A-6), respectively. 
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Figure Captions 

Fig. 1 The corrections of the initial gluons to the basic amphtude of the DGLAP 
equation and they lead to (b) BFKL-, (c) modified DGLAP- and (d) modified BFKL 
equations, respectively. 

Fig. 2 Application of the W — W approximation. The dark circles denote perturbative 
QCD interaction with the correlation of the initial partons. 

Fig. 3 A sub- vertex in a complex Feynman diagram. 

Fig. 4 The parton distributions on the impact space of a nucleon: (a) the dilute parton 
system, (b) some of partons form the cold spots (dashed circles) and (c) the saturation, 
where the long- and short-distance correlations among partons coexist in a nucleon. 

Fig. 5 The QCD correlations between two initial gluons. 

Fig. 6 The TOPT-diagrams consisted by the elemental amplitudes in Fig. lb on the 
impact space. These diagrams lead to the BFKL equation. The dashed hues are the time 
ordered lines in the TOPT and the dark box represent a non-local measurement of the 
unintegrated distribution. 

Fig. 7 The illustration of the right-hand side (real part) of (a) DGLAP equation, (b) 
BFKL equation, (c) modified DGLAP equation and (d) modified BFKL equation, where 
un-connecting gluons imply all possible connections inside the cold spot like in Fig. 6. 
The last sub-diagrams are added probe, which contribute S{x — 2 — xb) or Eq. (2.10) in 
(a) and (c), or in (b) and (d), respectively. 

Fig. 8 The virtual diagrams corresponding to Fig. 6. 

Fig. 9 A cutting diagram originating from Fig. Id. 

Fig. 10 The TOPT-diagrams consisted by the elemental amphtudes in Fig. Id in the 
impact space. These diagrams lead to the modified BFKL equation. 

Fig. 11 Two cutting diagrams contributing to the modified DGLAP equation. 



Fig. 12 The factorization of Fig. 10, which leads to an approximate form of the 
modified BFKL equation. 

Fig. 13 The unintegrated gluon distribution J[x^^,x) at = 5GeV~^ as a solution 
of the modified BFKL equation (4.26) taking the approximation (4.27) (sohd curve), a 
corresponding solution of the BFKL equation (dashed curve) and a solution of Eq. (4.26) 
but without antishadowing (point curve) . The broken-point line is an imaginal saturation 
limit. 

Fig. 14 A schematic kinematic region of four evolution equations in this work. 

Fig. 15 One-step evolution including the 4 — > 2 recombination amplitude, which 
is regarded as a graphic description of the nonlinear terms in the Balitsky-Kovchegov 
equation. 

Fig. 16 Schematic cut diagrams for the BFKL equation (a-d) and corresponding ladder 
diagrams (e-h). 

Fig. 17 A typical ladder diagram in the traditional BFKL theory, where the vertical 
(dashed) lines are reggeized gluons. 

Fig. 18 A comparison of our partonic picture with the color dipole approaches. 

Fig. 19 (a) A cut diagrams of the modified DGLAP equation and (b) corresponding 
ladder diagram. 

Fig. 20 The helicity amplitudes contributing to the spin-independent BFKL kernel, 
where the processes in (al)-(a4) should be inhibited at the W — W approximation. 

Fig. 21 Complete TOPT diagrams containing the probe for the DGLAP equation: (a) 
real diagram and (b) virtual diagrams. 



